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Abstract. Tlie basic tiieory of pliotoemission, inverse pliotoemission, Auger-electron 
and appearance-potential spectroscopy is developed within a unified framework starting 
from Fermi's golden rule. The spin-resolved and temperature-dependent appearance- 
O ' potential spectroscopy of band- ferromagnetic transition metals is studied in detail. It is 

5^ , shown that the consideration of electron correlations and orbitally resolved transition- 

"^ ■ matrix elements is essential for a quantitative agreement with experiments for Ni. 
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1 Basic electron spectroscopies 



^ I Electron spectroscopy ||l|,plpPP,pl is one of the fundamental experimental tech- 

Q ■ niques to investigate the electronic structure of transition metals. For a deeper 

O I understanding of the electronic properties, a meaningful interpretation of the 

measured spectra is necessary which is based on a reliable theory of spectro- 
scopies. This is particularly important for the study of band ferromagnetism 
^ ■ which is caused by a strong Coulomb interaction among the valence electrons. 

t^ I Due to the presence of strong electron correlations, simple explanations of spec- 

tral features within an independent-particle model may fail. Here we try to clarify 
what physical quantity is really measured and what ingredients are needed for a 
(•~^ • theoretical approach to find a satisfactory agreement with the experimental data. 

We are concerned with the theory of four basic types of electron spectroscopy 
(see Fig. |l|): photoemission (PES), inverse photoemission (IPE), Auger-electron 
(AES) and appearance-potential spectroscopy (APS). 

a The spin-, angle- and energy-resolved (ultraviolet) valence-band photoemis- 

I sion (PES) more or less directly measures the occupied part of the band struc- 

i-rt ■ ture. Present theories of PES are mostly based on the so-called one-step model 

P^ I of photoemission |^,|3!||,^ which treats the initial excitation step, the transport 

O ■ of the photoelectron to the surface and the scattering at the surface barrier as a 

single quantum-mechanically coherent process. The one-step model is essentially 
based on the independent-particle approximation. The valence electrons move 
independently in an effective potential as obtained from band-structure calcula- 
tions within the local-density approximation (LDA) of density- functional theory 
C^ I (DFT) [|lO|,0. There are recent attempts lljjl^l for a reformulation of the one- 

step model to include a non-local, complex and energy-dependent self-energy 
which accounts for electron correlations. 

The one-step model also applies to inverse photoemission spectroscopy. IPE 
is complementary to PES and yields information on the unoccupied bands above 
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Fig. 1. Schematic picture of the initial and the final states for different electron 
spectroscopies. PES: photoemission spectroscopy (s = 1). IPE: inverse photoemis- 
sion (s = —1)- AES: Auger-electron spectroscopy (s = 2). APS: appearance-potential 
spectroscopy (s = —2). The grey region is the occupied part of the valence band 



the Fermi energy |9|. For example, IPE can be used to determine the dispersions 
of the Rydberg-like series of surface states in the 1/z image potential in front of 
the surface g]- 

The theory of CVV Auger-electron spectroscopy (AES) is more complicated 
compared to PES/IPE since there are two valence electrons participating in the 
transition. Due to core state involved additionally, AES is highly element spe- 
cific. This feature is frequently exploited for surface characterization. Contrary 
to k resolved (inverse) photoemission, the Auger transition is more or less lo- 
calized in real space. High-resolution AES may thus yield valuable information 
on the local valence density of states (DOS). In the most simple theoretical 
approach suggested by Lander in 1953 |14| , the Auger spectrum is given by 
the self-convolution of the occupied part of the DOS. Modern theories of AES 
also include the effect of transition-matrix elements ||l5|| . Similar as the one-step 
model of PES and IPE, these approaches are based on the independent-electron 
approximation. 

Finally, the appearance-potential spectroscopy (APS) is complementary to 
AES. Within Lander's independent-electron model [|^ the AP line shape result- 
ing from CVV transitions is given by the self-convolution of the unoccupied part 
of the DOS. Matrix elements are included in refined independent-particle theo- 
ries [|l6| . Its comparatively simple experimental setup and its surface sensitivity 
qualifies APS to study surface magnetism, for example. For a ferromagnetic 
material, the spin dependence of the AP signal obtained by using a polarized 
electron beam gives an estimate of the surface magnetization as has been demon- 
strated for the transition metals Fe and Ni |17| , p^ . 

Common to all four spectroscopies are a number of fundamental concepts 
and approximation schemes used in the theoretical description. Therefore, it 
should be possible to develop the basic theory of PES, IPE, AES and APS 
within a unified framework. This is discussed in the following section y. As a 
result of this basic theory of electron spectroscopies it found that in each case 
the intensity is essentially given in terms of a characteristic Green function. 
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The actual evaluation of the Green function is more specific and depends 
on the respective spectroscopy, on the material and the physical question to be 
investigated. In section we will therefore pick up an example and give a more 
detailed discussion of APS for band ferro-magnets. 

A notable quality of APS is its direct sensitivity to electron correlations: In 
the AP transition, two valence electrons are added to the system at a certain lat- 
tice site. These final-state electrons mutually experience the strong intra-atomic 
Coulomb interaction U. As a consequence the two-particle (AP) excitation spec- 
trum is expected to show up a pronounced spectral-weight transfer or even a 
satellite with a characteristic energy of the order of U [|l^;Q . This is an impor- 
tant difference between the one- (PES, IPE) and the two-particle spectroscopies 
(AES, APS). To demonstrate this sensitivity to correlations in APS we consider 
the spin- and temperature-dependent AP spectra of ferromagnetic Nickel as a 
prototype of a strongly correlated and band-ferromagnetic material in section ]^ 
The discussion of the calculations and the comparison with experimental results 
will show what presently can be achieved in the theory of electron spectroscopies. 
Some conclusions are given in section o. 

2 Theory 

The basic theory starts from the Hamiltonian H which shall describe the elec- 
tronic structure of the system within a certain energy range around the Fermi 
energy. In addition, we consider a perturbation Vg that mediates the respective 
transition. Vg is assumed to be small and will be treated in lowest-order pertur- 
bation theory. The index s distinguishes the different spectroscopies (see Fig. |l|). 
s = ±1,±2 stands for the difference between the number of valence electrons 
before and after the transition. 

In the photocmission spectroscopy (PES, s = +1) a valence electron is excited 
into a high-energy scattering state by absorption of a photon. The photoclcctron 
escapes into the vacuum and is captured by a detector depending on its spin, 
energy and angles relative to the crystal surface. Within the framework of second 
quantization, the perturbation can be written as: 

V+i = Y^ M^^ aic^ + h.c. (PES) . (1) 

c-y annihilates a valence electron with quantum numbers 7, and aj, creates an 
electron in a high-energy scattering state a. The matrix element is calculated 
using the electric dipole approximation. This is well justified in the visible and 
ultraviolet spectral range. Neglecting the term quadratic in the field as usual 
and choosing the Coulomb gauge, one obtains: 

M„^ = («|Aop|7) (s = ±1) . (2) 

Here p is the momentum operator and Aq is the spatially constant vector poten- 
tial inside the crystal. It can be determined from classical macroscopic dielectric 
theory. Note that we use atomic units with e = ttIo = fi = 1. 
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In an invcrsc-photocniission experiment (IPE, s = —1) an electron beam 
characterized by quantum numbers a hits the crystal surface. An electron may 
be de-excited into an empty valence state 7 above the Fermi energy by emission 
of a photon. The photon yield is measured as a function of a (the polarization, 
the energy and the angles of the incident electron beam). We have: 

y_i - Y. ^7" ""W + i^-c- (IPE) . (3) 

Obviously, V-i = Vi. Both, PES and IPE, are induced by the same electron- 
photon interaction. 

In the case of the two-particle (s — ±2) spectroscopies AES and APS, the 
transitions are radiationless. Here the Coulomb interaction mediates between 
the initial and the final state. The initial state for AES (s = -f 2) is characterized 
by a hole in a core level (produced by a preceding core-electron excitation via 
x-ray absorption, for example). In the Auger transition the empty core state P is 
filled by an electron from a valence state 7. The energy difference is transferred 
to another valence electron in the state 7' which is excited into a high-energy 
scattering state a and detected. Thus, 

V+2 - J2 ^"/37'7 aiblc^cy + h.c. (AES) . (4) 

a/37'7 

The creator bt refers to the core state. 

Finally, an electron (a) approaching the crystal surface can excite a core elec- 
tron (/3) into a state of the unoccupied bands. Above a threshold energy both, the 
de-excited primary electron and the excited core electron occupy valence states 
(7, 7') above the Fermi energy after the transition. The appearance-potential 
spectroscopy (APS, s = —2) monitors the intensity of this transition as a func- 
tion of polarization, energy and momentum of the incoming electron beam by 
detecting the emitted x-rays or Auger electrons of the subsequent core-hole de- 
cay. The perturbation is: 

^-2= II ^^7%a/3 4'4V«+ll.C. (APS). (5) 

a/37'7 

Again, V-2 = V2. Since AES and APS are induced by the Coulomb interaction, 
the matrix element reads: 

M„^y, = «(a|(2)(/3|^^|y)(i)|^)(2) (s = ±2). (6) 

To calculate the cross sections for the different spectroscopies, time-depen- 
dent perturbation theory with respect to the perturbation Vs is applied. Using 
Fermi's golden rule, the transition probability per unit time for a transition from 
the initial state \Ei) to the final state \Ef) is given by: 

Ws = 27r\{Ef\Vs\E,}\^ 6{Ef - E, - io,) . (7) 
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\Ei) and \Ef) are eigenstatcs of the grand-canonical Hamiltonian H — fiN with 
eigenenergies Ei and Ef . fi is the chemical potential and A^ the particle- number 
operator. Furthermore, lUs = ±a; for s = ±1 where lu is the photon frequency. 
{us =Ofor s = ±2). 

We proceed by applying the so-called sudden approximation which reads: 

\Ef) = ai\E^) , Ef =e„ + En {s = +1, +2) , 

\E,) = ailEn) , E,^e„ + En {s - -1, -2) . (8) 

At this point we have neglected the interaction of the high-energy electron with 
the rest system. The latter is left in the n-th excited state of H with eigenenergy 
En- f-a is the one-particle energy of the high-energy scattering state. The electron 
and the rest system propagate independently in time but consistent with energy 
conservation. Generally, the sudden approximation is believed to hold well if Eq, 
is not too small. We can furthermore assume that 

a„|^,)-0 (s = +l,+2), a^\Ef)~Q (s = -1, -2) . (9) 

Hence, 

{Ef\Vs\E,) = (^„| [a„,I4]_ \E,) = (^„|Ti^)|S,) {s = +l,+2) , 

{Ef\V,\E,) - {Ef\ [1/ ,at.]_ \E^) ^ (i?„|ri«) V.) {s = -1, -2) , (10) 

where [..., ...]_ denotes the commutator and 

ri*) = K, K]_ - Y. ^"-z '^"'^^ ('^ = ±1) - 

7 

T^^) = K, K]_ = Y. ^"/37'7 b^c^c^' {s = ±2) (11) 

the transition operator. The transition probability per unit time now reads: 

ws - 2^\{Er,\Ti'^\Ei)\' 5{En + e^-E,- lUs) (s = +1, +2) , 

Ws^2Tr\{Ef\T^'^^\E^)\H{Ef-E,,^e^-Lu,) (s = -1, -2) . (12) 

To get the intensity /<, we have to average over the possible initial states. 
Initially, the system is assumed to be in thermal equilibrium: At the tempera- 
ture T {(3 — l/fceT) the system is found in the state \Ei) with the probability 
Wi — Z^^ exp{—(3Ei). Here Z — ^^ exp{—f3Ei) is the grand canonical partition 
function. We consider an experiment that determines all quantum numbers a 
necessary for a complete measurement (s = -1-1, +2) or preparation (s = —1, —2) 
of the state of the high-energy electron. Consequently, all indices have to be 
summed over except for a. Eventually, this yields the intensity: 

^-(") = Y E^^'^'^'l^^"!^"''!^*)!" '^(^" +e^-E^- LUs) {S = +L+2) , 

i,n 

^«(") = ^E^"''''"K^/I^"'^V«)P'5(£;/ -i?n -^a -C..) {S = -1,-2) . 



z 



(13) 
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These expressions can be written in a more compact form. For this purpose 
we define the Green function [pl| as: 



m,n 



which can also be written in the form Gaa'{E) = {{Ta ,T^, )) to show the 
dependence on the transition operator. With the help of the Dirac identity 

1/{E + iQ+) = V{l/E) - iTTS{E), 

h{a) = -2 j} Im G„„(S + ^0+) (s = +1, +2) , 

h{a) = -2 .„ , ^ ImGo,a{E + iQ+) (s = -l,-2). (15) 



2/3S + 1 



E=ea 



One recognizes the Fermi function f{E) = l/{exp{/3E) + 1) and the fact that 
both, the direct (s > 0) and the inverse spectroscopies (s < 0) are described by 
the same (Green) function. In fact, the following relations hold: 

I+iiE)^e-P^I^i{E), I+2{E)^e-P^I^2{E). (16) 

Eqs. (|lj) and (|l5|) show that the intensity is given by a weighted sum of 5 peaks 
at the excitation energies E — En — Em- In the thermodynamic limit the energy 
spectrum will be continuous in general and thus the dependence Is {E) is smooth. 

The weight factors {E„\T^, \Em) {Em\Ta' |£'„) distinguish between the different 
spectroscopies. 

The final equation (|l^) is the goal of our considerations so far. It relates 
the intensity to the Green function (nj) which is a central quantity of many- 
body theory. It can be (approximately) calculated by standard diagrammatic 
methods such a perturbation theory with respect to the interaction strength or 
by methods involving infinite re-summations of diagrams [ pl| . This is an essential 
advantage compared with a direct evaluation of Eq. (|l^) . The latter seems to be 
impossible since one would have to compute explicitly the eigenenergies En and 
eigenstates |£'„) of a system of interacting electrons. 



3 APS for band ferro- magnets 

In the following we will concentrate on the appearance-potential spectroscopy 
to give an example how calculations can be performed in practice. Beforehand, 
however, some preparations are necessary. 

Core-hole effects. A characteristic feature of APS is the formation of a core 
hole. The energy of core level e^ involved in the transition determines the shallow 
energy: Energy conservation requires that the energy loss of the primary electron 
must be equal to or larger than Ep — ec where Ep is the Fermi energy. Besides 
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this static effect there are also dynamic core-hole effects originating from the 
scattering of valence electrons at the local core-hole potential in the final state 
for APS p3,p3. The dynamic core-hole effects are usually neglected assuming 
the Coulomb correlation between valence and core electrons to be small and 
not to affect the AP line shape significantly. This is an approximation which is 
difficult to justify and which has to be checked in each case separately. 

It leads, however, to a substantial simplification of the problem. Similar to 
the sudden approximation discussed above, one can write: 

\E,) ^ alhl\E„)\ E, = e,, + e, + E'„, bp\Ef) = . (17) 

With essentially the same steps as above one gets: Zaps oc (1 — ,f{E)) lTiiG{E + 
iO"*") where now E — ea + Cc and where G is a two-particle Green function of the 
type G = {{ J2 Mc^c-y' ;^M*cjJ,c!,,)). Note that the core-electron creators are 
eliminated (cf. Eq. (pi)). The only dynamic degrees of freedom left are those of 
the valence electrons (7). 

Hamiltonian. Characteristic for the valence electronic structure of the band- 
ferromagnetic 3d transition metals are the strongly correlated 3d bands around 
the Fermi energy which hybridize with essentially uncorrelated 4s and 4p bands. 
The band structure derives from a set of localized one-particle basis states I7) 
with 7 specified as 7 = {i,L,a). Here \iLa) is taken to be a localized (atomic- 
like) orbital centered at the site z of a lattice with cubic symmetry, a =1, i is the 
spin index. L is the orbital index running over the five 3d orbitals, the 4s and 
the three 4p orbitals. We also introduce an index m which labels the different 
d orbitals, namely the three-fold degenerate t2g and the two-fold degenerate Cg 
orbitals. Using these notations the Hamiltonian H reads: 

ii'LL'a iaa' m\...m4^ 

(18) 
This is a multi-band Hubbard-type model including a strongly screened on- 
site Coulomb interaction among the d electrons. The hopping term oc t^f" = 
{iLa\H{U... = 0)\i'L'a) describes the "free" (non-interacting) band structure 
which can be obtained by Fourier transformation to k space tfj^ i— > t^-^ (k) and 
subsequent diagonalization i^^ (k) t-^ er(k). 

AP intensity. Having specified the one-particle basis and the Hamiltonian, we 
can write down the final expression for the AP intensity with all relevant depen- 
dencies made explicit: 

/.,.,(k||,i?)cxlm Y. ^■^ra(k|h^)x 

x{{c^L,..C.L2..;cl,,AL',J)E (Af^f ^', (k|| , i?))* . (19) 

The intensity depends on the quantum numbers f3 of the core hole formed, 
particularly on the spin CTc of the core state, and on the quantum numbers a of 
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the incoming electron, its energy E, its momentum parallel to the crystal surface 
kii and its spin (7^. 

The AP line shape essentially results from intra-atomic transitions. Conse- 
quently, those transition-matrix elements M that lead to off-site contributions 
to the intensity are neglected in (ttW). The "raw spectrum" as given by the imag- 
inary part of the Green function in Eq. ( [19[ ) is therefore isotropic. Any angular 
(k||) dependence is due to the angular dependence of the high-energy scattering 
state in the matrix element. 

The energy dependence of the intensity, i. e. the actual AP line shape, is 
mainly determined by the two-particle Green function and reflects the energy- 
dependent probability for two-particle excitations. On the contrary, for typical 
kinetic energies of the primary electron of the order of keV the change of the 
matrix elements due to the energy dependence of the high-energy scattering 
states is expected to be weak over a few eV. 

Spin dependence. The Coulomb interaction that induces the transition conserves 
the electron spin orientation. So the spin orientations of the incoming electron 
ai and of the core electron ac in the initial state determine the spin orientations 
of the two additional valence electrons in the final state. For an incoming elec- 
tron with spin orientation ai one can distinguish between a "singlet" transition, 
i. e. excitation of a core electron with ac — —ai, and a "triplet" transition with 
ac — ai. It is important to note that the ratio between singlet and triplet tran- 
sitions is regulated by the symmetry behavior of the transition-matrix elements 
under exchange of the orbital indices. Assume, for example, that the matrix el- 
ement is symmetric: Ml-^l^ = M^^Li. Then the transition operator for triplet 
transitions vanishes, T^^cn = J2liL2 ^^LiL2 CtLia^CiL^a, = for at = ac, since 
cf^g. = (Pauli principle) and since CiL^cr^CiL2a-i is antisymmetric with respect to 
Li <^ L2. Thus, a symmetric or even a constant matrix element (as is sometimes 
assumed for simplicity) completely excludes the triplet transitions. 

Consider a paramagnetic material for a moment. The argument above shows 
that a symmetric M^-^l^ would imply ac — —ai, i. e. a fully polarized primary 
electron beam leads to a full polarization of the core hole. Thus, any deviation 
from full core-hole polarization must be due to the antisymmetric part of the 
matrix element. 

Below we will consider a ferromagnetic material and a situation where the 
spin state of the final core hole is not detected. Then, the intensities have to 
be summed incoherently: la-^ = la^ai + I-a^ai- For a ferromagnetic material 
one expects the intensity to be still dependent on the spin orientation of the 
primary electrons: /| 7^ /^ (while /| — /| above the Curie temperature). Let us 
assume again symmetric behavior: M^-^l^ — Ml^l-^. It has been argued above 
that this implies la-ao-i = for ai ~ ac. Furthermore, it is easy to see that 
T|| = — Tji which implies /|| — /||. Hence, /| = /j- ^^ conclusion, any spin 
asymmetry in the intensity is due to a non-vanishing antisymmetric part of the 
matrix elements. The AP intensity asymmetry is much more determined by the 
symmetry properties of the matrix elements with respect to their orbital indices 
as compared to their spin dependence. 
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Direct and indirect correlations. Any diagrammatic approach gives the two- 
particle Green function {{cc;c'^c^)) as a (comphcatcd) functional J-" of the one- 
particle Green function ((c;ct)). One can thus distinguish between direct and 
indirect correlations [e3. As has been shown above, the one-particle Green func- 
tion corresponds to the (inverse) photoemission spectrum. Indirect correlations 
in APS are those which originate from the renormalization of the free one-particle 
spectrum by the interaction term in (jTsI). The direct correlations, on the other 
hand, are represented by the concrete form of the functional T. Essentially, the 
direct correlations originate from the direct Coulomb interaction of the two ad- 
ditional valence electrons in the final state. Since both electrons are created at 
the same site, they are affected by the strong intra-atomic interaction. This may 
give rise to correlation-induced satellites in the AP spectrum as is demonstrated 
by the so-called Cini-Sawatzky theory ||l^,0 . 

When neglecting the direct correlations, the functional J-' becomes a mere 
self-convolution of the one-particle Green functions. If furthermore the transi- 
tion-matrix elements are taken to be constant, the AP spectrum is simply given 
by the self-convolution of the unoccupied part of the one-particle density of 
states oc Im((c; c^)). This is the so-called Lander model ||lj] which is frequently 
employed for a rough interpretation of the spectra. 

Within the framework of the self-convolution (Lander) model the two final- 
state electrons propagate independently. As a consequence one finds that only 
the direct term with Li = L{ and L2 = -^2 and the exchange term with Li = L'2 
and L2 = L'l contribute to the sum over the orbital indices. Generally, however, 
the Green function in Eq. (^^ depends on four orbital indices. This implies that 
the usual characterization of the final state with two quantum numbers {d-d, 
s-d, etc.) is no longer valid if the direct correlations are included. The orbital 
character may change by electron scattering. 

4 Appearance-potential spectra of Nickel 

The significance of electron-correlation effects in APS shall be elucidated in a 
more detailed discussion below. For this purpose we concentrate on ferromagnetic 
Nickel as a prototypical 3d band-ferromagnet and compare the results of the 
theoretical approach with experimental data. 

Experiments. Experimental results are available for a Ni(llO) single-crystal sur- 
face with in- plane magnetization [p4[|25|] . In the setup a spin-polarized electron 
beam is used for excitation which is emitted from a GaAs source. To correct for 
the incomplete polarization of the electrons {P w 30%), all data are rescaled to a 
100% hypothetical beam polarization. The spin effect is maximized by alignment 
of the electron polarization and the sample magnetization vector. The core-hole 
decay is detected via soft-X-ray emission (SXAPS). To separate the signal from 
the otherwise overwhelming background, modulation of the sample potential by 
a peak-to-peak voltage of 2V together with lock-in technique is employed. Details 
of the experimental setup can be found in Rcfs. |24,nN]. 
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Fig. shows the measured differential AP intensity as a function of the 
energy of primary electrons with polarization parallel (minority, J,) or antiparallel 
(majority, t) to the target magnetization. The displayed energy range covers the 
emission from the ini transition (2^3/2 core state). The Lu (2^1/2) emission 
would be seen at higher energies shifted by the 2p spin-orbit splitting of 17.2 eV. 

For T/Tq — 0.16 (Tc « 630 K) the system is close to ferromagnetic satu- 
ration. The AP spectrum shows a strongly spin-asymmetric intensity ratio as 
well as a spin splitting of the main peak a,t E = 852.3 eV (indicated by the 
dotted line). Since Ni is a strong ferromagnet, there are only few unoccupied 
d states available in the majority spin channel, and thus /| < I^ holds for 
the (non-differential) intensities. This is the dominant spin effect. The inten- 
sity asymmetry in the main peak gradually diminishes with increasing T and 
vanishes at Tq- 

The main peak is related to the high DOS at the Fermi energy. Within an 
independent-electron picture, one can thus characterize the main peak as orig- 
inating from transitions with two final-state electrons of d-d character mainly. 
This is corroborated by calculations based on DFT-LDA ||l^,|l6[. Additional 
small s-d contributions are present in the secondary peak at i? w 859 eV as has 
been concluded from the analysis of the transition-matrix elements. The sec- 
ondary peak has been identified as resulting from a DOS discontinuity deriving 
from the L7 critical point in the Brillouin zone [ [L8| . No temperature dependence 
and spin asymmetry is detectable here. 

Hamiltonian. To study the significance of electron correlations, we consider a 
nine-band Hubbard-type model H — Hq + Hi including correlated 3d and uncor- 
related 4s/4p orbitals as given by Eq. ([18|). The first ("free" ) term Hq is obtained 
from a Slater-Koster fit to the paramagnetic LDA band structure for Ni Pqj . 
Opposed to PES/IPE, this comparatively simple tight-binding parameterization 
appears to be sufficient in the case of APS since the two-particle spectrum does 
not crucially depend on the details of the one-particle DOS. 

The on-site interaction among the 3d electrons is described by the second 
term Hi. Exploiting atomic symmetries the Coulomb-interaction parameters 
which depend on the four orbital indices mi, . . . , m^ can essentially be expressed 
in terms of two independent parameters U and J. The numerical values for the 
direct and exchange interaction, U = 2.47 eV and J = 0.5 eV, are taken from 
Ref. pfl] where they have been fitted to the ground-state magnetic properties of 
Ni. Residual interactions involving de-localized s and p states are assumed to be 
sufficiently accounted for by the LDA. Finally, a double-counting correction is 
applied to H since partially the interaction Hi is already included in the LDA 
Hamiltonian Hq (for details see Ref. [p6|). 

Green functions. The two-particle Green function in Eq. ( |l9|) is approximately 
calculated by using standard diagrammatic techniques. Because of the low den- 
sity of 3d holes in the case of Ni, it appears to be reasonable to employ the 
so-called ladder approximation ||27| which extrapolates from the exact (Cini- 
Sawatzky) solution for the limit of the completely filled valence band [H 
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Fig. 2. (from Ref. |G3]) Spin- resolved Lm Ni AP spectrum for different reduced tem- 
peratures T/Tc- Data points: measured differential intensity dl /dE as a function of the 
primary energy. For better comparison with theory the same data are shown twice (left 
and right panel) . Lines, left: theory with direct and indirect correlations included (lad- 
der approximation). Lines, right: indirect correlations included only (self-convolution) 



For finite hole densities the ladder approximation gives the two-particle as a 
functional J- of the one-particle Green function (direct correlations). 

The one-particle Green function, which describes the indirect correlations, 
is calculated self-consistently within second-order perturbation theory (SOFT) 
around the Hartree-Fock solution |2g] . For a moderate U and a low hole density, 
a perturbational approach can be justified p8[ . A re-summation of higher-order 
diagrams is important to describe bound states ( "Ni 6 eV satellite" ) [g^ which, 
however, are relevant for AES only. Since spin-wave excitations are neglected in 
the approach, the calculated Curie temperature Tc = 1655 K is about a factor 
2.6 too high. Using reduced temperatures T/Tq, however, the temperature trend 
of the magnetization is well reproduced p6[ . 



(20) 



Matrix elements. The transition- matrix elements in Eq. (|19[), 

M^J^;(k||,^) = (2p,a,| (k||Sa,| r^; \iL,a,) \1L2a,) , 

are calculated by assuming the transition to be intra-atomic as usual [pl jl^jl^ . 
The different wave functions as well as the Coulomb operator l/ri2 are expanded 
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into spherical harmonics, the angular integrations are performed analytically, 
and the numerical radial integrations are cut at the Wigner-Seitz radius. 

Surface effects enter the theory via the high-energy scattering state \k^^Eai). 
It is calculated as a conventional LEED state with k|| = to describe the nor- 
mally incident electron beam in the experimental setup. The (paramagnetic) 
LDA potential for Ni is determined by a self-consistent tight-binding linear 
mufRn-tin orbitals (LMTO) calculation p^. The 3d, 4s, and 4p valence orbitals 
\iLa) are taken to be the muffin-tin orbitals. The four-fold degenerate 2p3/2 
core state is obtained from the LDA core potential by solving the radial Dirac 
equation numerically. Its (relativistically) large component is decomposed into 
a (coherent) sum of PauH spinors |2p, CTc) with CTc =t, i- 

Results. The solid lines in Fig. g (left) show the spectra as calculated from 



Eq. (19) using the ladder approximation for the two-particle Green function. To 
account for apparatus broadening, the results have been folded with a Gaussian 
of width a = 0.6 eV (see Ref. [Q). The calculated data are shifted by 852.3 eV 
such that onset of the un-broadened spectrum for T/Tq — 0.16 coincides with the 
maximum of im emission in the experiment (dotted line). Fig. y and also a more 
detailed inspection show that the secondary peak at E' « 859 eV is not affected 
by correlations at all. This is consistent with observed temperature independence 
of the peak and with the fact that the DOS has mainly s-p character at the 
discontinuity deriving from the Ly critical point. The maximum of the secondary 
peak is therefore used as a reference to normalize the measured spectra for each 
temperature. 

What are the signatures of electron correlations? The indirect correlations 
manifest themselves as a renormalization of the one-particle DOS. Here, they are 
responsible in first place for the correct temperature dependence of the intensity 
asymmetry of the main peak in the AP spectrum. Details are discussed in Refs. 
[E6|23] . The the direct interaction between the two additional final-state electrons 
and thus the direct correlations are much more important for APS. To estimate 
this effect, Fig. (right) also displays the results of the self-convolution model 
for comparison (still including matrix elements as well as the fully interacting 
one-particle DOS). 

Looking at the results of the ladder approximation, the overall agreement 
with the measurements is rather satisfying. Except for the lowest temperature 
the intensity, the spin splitting and the spin asymmetry of the main peak are 
well reproduced and, consistent with the experiment, a negligibly small inten- 
sity asymmetry and spin splitting is predicted for the secondary peak. However, 
switching off the direct correlations (Fig. ||, right), results in a strong overesti- 
mation of the main peak structure. 

A plausible qualitative explanation of this pronounced correlation effect can 
given within the Gini-Sawatzky theory: For low hole density the main effect of 
the direct correlations is known to transfer spectral weight to lower energies 
inaccessible to APS. This weight shows up again in the (complementary) Auger 
spectrum (recall that APS and AES are described by the same Green function). 
Hypothctically, for U ^^ oo all weight would be taken by a satellite split off at the 
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Fig. 3. (from Ref. |||l) Ni AP spectrum for T = 0. (a) full theory, (b) as (a) but 
matrix elements taken to be constant (see text) 



lower boundary of the Auger spectrum ||l9|,^ . A considerable weight transfer is 
in fact seenin AES |,|l|. 

Fig. shows the effect of the transition-matrix elements. Their importance 
for a quantitative understanding of spin-resolved APS can be demonstrated by 
setting M'^"'^' (k||,i?) = ±1 = const for Li > L2 or Li < L2, respectively, 
(see below) and comparing with the results of the full theory. Their energy 
dependence (via the energy dependence of |k||i?CTi)) is weak over a few eV at 
energies of the order of keV and cannot explain the difference between (a) and 
(b) in Fig. ^. The main difference is rather a consequence of the fact that the 
radial 2p core wave function has a stronger overlap with the (more localized) 
3d as compared to the (more de-localized) 4s/4p radial wave functions. This 
implies a suppression of the s-p contributions to the orbital sum in Eq. (|l9|). 
The features above E = 860 eV originate from additional discontinuities of the 
s-p-like DOS (as for the peak at £; « 859 eV). 

For T < Tq the spin asymmetry of the spectrum is mainly due to the spin 
dependence of the Green function in Eq. ( [l9| ) . If the calculation of the matrix 
elements (^0|) starts from the spin- polarized L(S)DA potential, an additional spin 
asymmetry is observed resulting from the spin dependence of the states in Eq. 
(po|). This, however, is small and has practically no influence on the results. 

On the other hand, Fig. y shows a strong suppression of the intensity asym- 
metry at high energies when taking matrix elements into account. This effect 
is controlled by the symmetry of the matrix Ml^l^ = M^''^'^(k||,_E). In the 
antisymmetric case, Ml-^l^ = — M^^Li {Li ^ L2), there is a maximum spin 
asymmetry (Fig. |^) while, even for a ferromagnet, there is no spin asymmetry 
at all for the symmetric case, as discussed in section ||. The results of the full 
calculation along Eq. (|20|) are neither fully symmetric nor antisymmetric with 
respect to Li, L2- 
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5 Conclusions 

It has been shown that the basic theory of different electron spectroscopies, 
PES, IPE, AES and APS, can be developed within a unified framework. Starting 
from Fermi's golden rule and employing the sudden approximation, the intensity 
can be expressed in terms of a one- or two-particle Green function including 
electron-photon or Coulomb matrix elements, respectively. Citing APS from the 
typical band ferromagnet Ni as an example, it could be demonstrated that a 
quantitative agreement with spin-resolved and temperature-dependent measured 
spectra can be achieved when the basic spectroscopic formula is taken seriously 
and is evaluated by modern techniques of solid-state theory. 

In particular, it has been demonstrated that there are pronounced (direct) 
correlation effects in the AP line shape of Ni. While s-p derived features at higher 
energies appear to be sensitive to the geometrical structure only, the main peak 
is strongly affected by the direct interaction between the two additional final- 
state electrons. Consistent with the Cini-Sawatzky model and consistent with the 
well-known Ni Auger spectrum, there is a considerable spectral-weight transfer 
to energies below the threshold. For Co and Fe one can even expect stronger 
effects of d-d correlations on the AP line shape since the d-hole density is larger 
than in Ni. Simple self-convolution models neglecting the direct correlations must 
be questioned seriously. 

It has also been demonstrated that a mere computation of the (two-particle) 
Green function is insufficient to describe spin-resolved APS. The spin asymmetry 
of the AP signal is found to be mainly determined by the orbital-dependent 
transition-matrix elements and their transformation behavior under exchange of 
the orbital quantum numbers. 

Considering the temperature dependence of the spectra and the magnetic 
order resulting from strong (indirect) correlations in addition, one can state that 
the AP line shape of a typical ferromagnetic Zd transition metal is determined 
by a rather complex interplay of different factors. 

Despite the fact that a reasonable understanding of APS from Ni has been 
achieved, there is much work to be done in the future: An open question con- 
cerns the importance of core-hole effects in APS, for example. For the present 
case there has been no need to consider scattering at the core-hole potential in 
the final state. This may likely be different for systems with a smaller id oc- 
cupancy. Furthermore, one must recognize that even the determination of the 
valence-band Green function is a central problem of many-body theory. Recently 
much progress has been achieved to deal with single-band models [^ ; for a quan- 
titative interpretation of electron spectra from transition metals, however, one 
needs a theory that realistically includes orbital degeneracy and sp-d hybridiza- 
tion from the very beginning. Compared to perturbation theory or lowest-order 
re-summation of diagrams, as employed here, improvements are conceivable and 
necessary albeit not easily performed. 
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